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Abstract 

We show how the operator product expansion evaluated in the approximation of ignor- 
ing gluons leads to the covariant formulation of the quark parton model. We discuss the 
connection with other formulations and show how the free quark model prediction, §2 = 0, 
changes smoothly into the Wandzura-Wilczek (WW) relation for quark masses small rel- 
ative to the nucleon mass. Previous contradictory parton model predictions are shown 
to follow from an inconsistent treatment of the mass shell conditions. The description is 
extended to include quark mass corrections. 



1 Introduction 



A recent preliminary measurement of the proton structure function (72 IQ in polarised deep 
inelastic scattering suggests that the integral form of Wandzura-Wilczekp], || relation 

92{x) = -gi{x) + / —gi{x') (1) 

is consistent with the new data — at least for x > 0.2. Here gi is the structure function already 
measured with good accuracy from previous experiments [H, ^, |^ with longitudinally polarised 
protons. Polarising the proton in the transverse direction gives a measurement of gi + g2 and 
hence an estimate for g2 can be extracted from both sets of data. 

In the operator product expansion (OPE) analysis of deep inelastic polarised scattering 
both twist-2 and twist-3 operators contribute to g2 while only the former contribute to gi. 
If one assumes the twist-3 operator contributions are negligible then the moments of (72 are 
related to those of gi and, assuming the Mellin transform can be performed]^, the relation of 
eq([^) appliesip. The motivation for ignoring the twist-3 terms follows from the observation]^ 
that, in the limit of zero quark mass, all twist-3 operators involve gluons. Thus in the quark 
model, with massless quarks and where gluons are neglected, one may expect the WW relations 
to apply because the twist-3 operator matrix elements vanish. In a previous paperf^ (JRR) 
we showed that a consistent covariant formulation of the quark parton model, in which the 
massless on-shell quarks have a non-zero transverse momentum kx, leads to precisely the WW 
relation eq(0). However this formulation has been questioned f^, because it is apparently 

in contradiction with other formulations of the parton model. 

Because of this and motivated also by the encouraging results of the preliminary SLAG 
data[Q which lends some support to the quark model approximation, we reconsider the deter- 
mination of the polarised structure functions in the quark model approximation i.e. neglecting 
the gluonic component of the nucleon but allowing for valence and sea quarks. The connec- 
tion with the OPE is discussed in detail and we show that the assumption of negligible gluon 
component in the contributing operators is equivalent to a covariant parton model description 
of the process in which the quarks, which may be massive, are allowed to have a completely 
general transverse momentum distribution. We find the covariant quark parton model provides 
^This may not be allowable due to singular (Regge) behaviour at small x in which case the moment relations 
only may apply Q. In what follows we ignore this potential problem. 



an entirely consistent picture of polarised scattering in which the form of g2 is fully determined 
once gi is measured. The origin of previous contradictions between quark model relations is 
shown to be due to incorrect treatment of the equations of motion in the associated OPE. 

We go on to consider the inclusion of quark masses in the parton model. In this case the 
parton polarisation vector picks up a component in addition to the longitudinal one and the 
WW relation is no longer satisfied. There is an additional polarised parton density associated 
with this component and restoring a connection between gi and g2 rests upon relating the 
parton densities associated with the two components. By considering plausible models for the 
source of the parton polarisation, we show how such a relation can be made, leading once again 
to a prediction of g2 in terms of the measured gi. 



2 Covariant parton model for on-shell massless partons 

We start with a brief review of the covariant parton model formulation of polarised scattering 
(in the massless limit) which we shall show is equivalent to the OPE in the quark model ap- 
proximation. The anti-symmetric hadronic tensor is written in terms of the two spin dependent 
structure functions gi and g2, 

K{q,P, S) = ie,,,a ^ [9i{x)S'^ + 92{x){S'^ - p"^)] (2) 

where Af, p, S arc the mass, momentum and polarisation vector of the proton; q is the momen- 
tum of the virtual photon and u = p.q/M. The covariant parton model expresses W^^ in terms 
of a convolution over the struck parton's momentum. 

Wl^f{q,p, S) = J2[ d'k Mp, k, S) w^M, k) 5[{k + qf - m'] (3) 

h ■' 

Here k,m,h are the parton's momentum, mass and (in the massless limit) helicity {h = ±1) 
and 

w^^ = ihme^^p^qPs" (4) 

where s"^ is the parton spin vector which as — > gives — k„. Summing over parton 
helicities thus gives the combination of parton densities f\ — f_ = Af{p,k,S) and since we 
insist on covariance, this can be written as Af{p.k,k.S,k^). As this combination has to be 
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linear in the proton polarisation (spacelike), covariance demands that A/ oc k.S and so we 
can write 



Afip.k, k.S, = J-t^f{p.k, e) (5) 
Substituting into eq(|^), comparing with eq(@) and using d^kdlik + g)^ — rn?] = {j^M / Au)dk'^dy 



k 

where y = 2^ = x ^ — ^ gives (for details see JRR) 



9i{^) = — Z — / "-y (2a; - y) h[y) (6) 



g2{x) = [' dy {2y - 3x) h{y) (7) 



where h{y) = J dk'^f{y, k"^). Crucial to deriving expressions (P,|3) and to the consistency of the 
parton model is the retention of the parton transverse momentum k^. In fact, from (|^J^) we see 
that the combination gi{x) + g2{x) (which is the relevant quantity for a transversely polarised 
proton) is proportional to {y — x) which is simply k^/xM'^. There is just one function h{y) 
(the relevant charge weighted combination of parton distributions) which determines both gi 
and g2 which is the key to obtaining a simple relation between the two. Differentiating the sum 
gi + g2 given by (|^J^) gives the relation eq(|l]). This can be alternatively expressed in terms of 
moments 

ri rn — ^ 1 

-0 (8) 



1 ^ 



■gi{x) + g2{x) 



n 

the Burkhardt-CottinghamfT^ (BC) sum rule being the n = 1 version 

nl 

dx g2{x) = (9) 

In section 4, we shall consider the corrections arising from allowing the quarks to have a 
non-zero mass. These corrections lead to a violation of the WW sum rules eq(|^) but making 
plausible assumptions for the origin of the parton polarisation, g2 can still be related to gi. The 
BC sum rule eq(^ is preserved exactly however. For a particular model the magnitude of the 
violation of the WW sum rules can provide a possible phenomenological estimate for the light 
quark mass. 



3 Consistency of the covariant quark parton model within 
the OPE 

At a more formal level the polarised deep inelastic scattering may be analysed using the OPE. 
Of course any viable model must be consistent with the OPE and here we show that the 
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covariant parton model is equivalent to the OPE approach in the assumption of neglecting the 
gluon component of the nucleon. Let us start with the operator product expansion description 
of polarised scattering 



n=0,2,... I i 

+ T.F3%^" - ^exo,^^')OZr''"{^.')] e,,...^ (10) 

i 

where the twist-2 operators are given by 

Q^^n.-f^n ^ i'^Si^{0)^'''y^D^'\..D^'-Xk^{0) - (Traces) (11) 

with Afc=i,.8 the SU (3) Gell-Mann matrices and -D^ = + igA^ is the covariant derivative. Si 
denotes symmeterisation in the cr/ii.../i„ indices. In addition there is a flavour singlet operator 
corresponding to the replacement of by the unit operator together with a flavour singlet 
operator involving gluon fields only that we do not give here. The twist-3 operators are given 

by 

_ r^A7/^(0)7V^'''-^''"Afc^(0) - [traces) (12) 

where S denotes symmeterisation of the /ii.../i„ indices and A anti-symmeterisation of the 
(TyUj indices. These operators have a (suppressed) renormalisation scale (yU^) dependence. The 
operator matrix elements are given by 

< P, S\02,k\p, S >= ■^^^^Si^S^Pi.^.-.p^^ - traces) 

< p, S\03^k\p, S >= ^2 1 SA{SaP^^^■■■P^^„ - traces) (13) 
Using this in eq(|IUp gives 

/ dx x" gi{x, Q"^) = -an, n = 0,2, 4... 
JO 4 

n 

/ rfxx"c/2(a;,Q^) = -— — -(rf„-a„), n = 2,4,... (14) 

In the case the twist-2 operators dominate, (i„ = leaving just one unknown matrix element 
per moment to describe two structure functions. Thus one obtains the WW relations eq(|^) for 
g2. Here we wish to demonstrate that these relations follow in the quark model simply by 
analysing the implications of the OPE under the assumption of negligible gluon content of the 
nucleon. Of course there is no guarantee that this assumption is realistic but it is of interest 
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to determine its phenomenological implications so that experiment may determine whether a 
pure quark model picture ever makes sense. Further, having a quark parton model (albeit 
equivalent to the OPE) is useful in guiding one's intuition when discussing the long-distance 
effects as we discuss below. One thing to note is that, due to QCD corrections, the pure quark 
model approximation can make sense only for one value of because the QCD corrections 
that give perturbatively reliable information about the evolution in require that the gluonic 
component must become significant at high even if the quark picture is adequate at low Q^. 
However if the pure quark model is a good approximation at one scale then one may conclude 
that the leading twist operators dominate at all higher scales, the gluonic component generated 
during evolution being entirely leading twist in this case. 




(a) (b) 

Figure 1: a) Form of the operator matrix elements in the absence of gluons. b) Form of the structure 
function using the OPE and neglecting the gluons both in the coefficient functions and the operator 
matrix elements. 

Given this motivation we now consider the form of the quark (and anti-quark) operator 
matrix elements corresponding to Fig. la. The contributing operators in this approximation are 
simply determined using the quark model to determine the coefficients -F2 3 of eq(|l3). Using 
this, the sum over operators with the operator matrix elements determined in the approxima- 
tion of neglecting gluons is equivalent to evaluating the diagram in Fig. lb. The latter requires 
knowlege of the distribution fs{p,S,k,s) = fs{p'^,k'^,p.k,k.S,p.s) of quarks with polarisation 
s within a proton. In the next section we will specify these polarisation states precisely but 
here we wish to continue our discussion within the framework of the OPE. The quark distribu- 
tions fsip^i k'^,p-k, k.S,p.s) in fact provide a convenient summary of the effects of the operator 
matrix elements and provide the easiest way of determining the implications of the valence 
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approximation for the general OPE. Let us note the following facts concerning the OPE: 

Operators related by equations of motion are not independent. In fact Shuryak and Vainsteinj^ 
used the equations of motion to show that, in the absence of quark masses, the twist-3 operators 
introduced above are equivalent to operators involving gluons only. Their general result (for 
massless quarks) is 

^ 1=0 
1=0 

where Gap = leapxaG^" . 

Thus the twist-3 operator matrix elements will vanish in the valence quark approximation 
up to mass effects. We shall discuss the massive case shortly; in any case we may use the 
equations of motion to set fc^ = m^, the quark mass squared. 

The second point is that the normalisation scale dependence of the operators is cancelled 
in the prediction for a physical quantity by the scale dependence of the coefficient functions 
^2^3. Thus the evaluation of Fig. lb which corresponds to the full observable cross section has 
no scale dependence. 

Putting these points together we see that the quark distributions in fact are functions of just 
two variables fsip"^, k'^,p-k, k.S,p.s) = fs{p.k, k.S) since p"^ and k"^ are given by the nucleon mass 
squared and the quark mass squared respectively and, in the relativistic quark limit, p.s = p.k. 

Let us discuss the implications for cases of increasing complexity. First consider the trivial 
case of scattering off a free on-shell quark moving collinearly with the target. In this case the 
quark distribution is known (trivial) and the scattering amplitude is simply given by the jq 
Born diagram. Some care must be exercised in the massless limit. Since g2 is to be determined 
via eq(0) (in this case p'^ and S'^ are the quark momentum and spin) we must use a massive 
quark with mass m otherwise the factor multiplying g2{x) vanishes. The massless quark limit 
is then easily obtained by taking the limit m — 0. The Born diagram gives an amplitude 
proportional to 

^me^,p^q's''6{{k + qy) (16) 
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Comparison with eq(0) immediately shows that (?2 = 0. How does this relate to the OPE 
analysis? The operators contributing simply correspond to those in the structure 

j^^^^r^f^M+4+^)i.i' (17) 

The term iplfj,{<i)lu'4' = 'i^fiauaQ^' '4^151'^'^ when combined with the factors found expanding (p + 
= ^i{—'2,p.q/q^y/ and Fourier transforming to change the to a derivative acting on the 
quark field, immediately leads to the combination of operators Q'^^^^■■■^^^ _|_ (^^Mi --/^n The reason 
this combination arises is simply because there is no symmeterisation (or antisymmeterisation) 
in the Born term between the a index and the indices associated with the {p.qY terms. As 
discussed above O2 alone gives rise to the WW relations and non- vanishing (72- Thus the 
vanishing of (72 must come about through a cancellation of the contributions of O2 and O3. At 
first sight this seems impossible because of the Shuryak Vainstein relation eq(15) which suggests 
O3 vanishes in the absence of gluons. However here this is not true as a careful computation 
of the mass effects reveals. To illustrate the point consider the simplest operators 

By writing 9^ = {7^,(^}/2 we may re-express these as 

O3 = imSAi^'y^-f^'yf.^df.^ij (19) 

We now take matrix elements of these operators between on-shell partons with polarisation 
vector Sfj_. In momentum space we have 



(18) 



< p, S\02\p, S > cc Tr{(^ + m)(l + 75'^)7a75^7w)PM2 

= Am^S^p^^ + S^^p^)p^2 

< P, S\Os\p, S > cc Tr{(^ + m)(l + 75'^)7a757m)PM2 

= Am{SaPf,, - S^,p^)p^^ (20) 

Thus we see explicitly that for the combination {O2 + O3) there is a cancellation between O2 
and O3 of the term oc Po- which is associated with g2. This cancellation persists in the massless 
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limit. The reason is that although the O3 (and O2) matrix element is formally of 0{m) this is 
cancelled by the 1/m behaviour of the polarisation vector ~ p^/ m. Thus O3 does contribute 
even in the limit gluons are ignored. We have somewhat belaboured our discussion of this very 
simple model because this point has been the source of considerable confusion in the literature. 

So far our discussion has shown that (72 vanishes in the free quark model, a well known result. 
Of course a model of free collinear quarks is not realistic and to determine the implications 
of more reasonable models we consider now the case of Fig. lb with nontrivial fs- As we 
have stressed the calculation of the full diagram is completely equivalent to the use of the 
OPE with operator matrix elements determined in the valence approximation. Following the 
analysis presented above the result is given in eqs. (||) and (0) and satisfies the WW relation. To 
understand the origin of this relation we note that the intermediate (vertical) states of Fig. lb 
are on-shell quarks and the result of eqs.(^ and (^ was derived to leading order only in the 
quark masses. To this order the quark polarisation vector is ~ k^/m and the matrix element 
of O3 between the quark states vanishes simply because one cannot form an antisymmetric 
tensor from alone. This immediately implies that only O2 contributes giving rise to the 
WW relation. This result again is puzzling because there is no obvious limit in which the free 
quark result with vanishing g2 can be obtained i.e. why doesn't this argument apply too in 
the free quark case discussed above? The answer is that the determination of (72 requires the 
identification of the coefficient of the second tensor in eq(0). As we noted above this vanishes 
identically in the free quark limit (the approximation that 5*^ ~ p^/ m) and for this reason we 
had to compute g2 in the massive quark case and take the massless limit only at the very end. In 
the present case however the tensor structure of eq(|^) refers to the real nucleon and therefore 
doesn't vanish when using the approximation ~ k^/m. Thus the argument showing the 
vanishing of the O3 contribution applies in this case and we do indeed get the WW relation. 
Of course in the limit the nucleon mass is taken small and the nucleon polarisation satisfies 
Sfj. ~ P/i/M the argument that O3 vanishes fails and the WW relation would not apply. This is 
the limit that establishes the connection with the free quark model result for it leads to (72 = 
again. Thus we may see the WW result applies in the limit m/M and m/kT are small while 
the free quark model result applies when m ~ M, with m/kT small. In the next section we 
extend the calculation to include mass effects to approach the region where m/M and m/kx 
are not negligible. 

Before leaving the OPE analysis we should comment on the role of another twist-3 operator 
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not included in the above analysis. This is the operator 

Om = mSAi)-f„-f^d^^...d^^ij (21) 



From eq ([l7|) we see this operator arises through the term proportional to m. However it is 
easy to check in the free quark model case that it gives rise to a gauge variant term which is 
cancelled by the term proportional to ^ in eq(|l^which involves the operators 02,3 and which 
have a matrix element proportional to m. 

Our analysis of the OPE in the approximation of ignoring gluons has shown that it is 
equivalent to a covariant formulation of the quark parton model. As discussed in this 
leads to a completely consistent description of both polarised and unpolarised scattering. The 
analysis also shows why other parton model formulations have led to contradictory results 
because the OPE requires the use of equations of motion to relate the various operators. This 
is in conflict with a parton model formulation in which the parton momentum is a fraction, x, 
of the nucleon momentum because the parton mass will then be xM, an x dependent quantity 
not consistent with the equations of motion of the underlying quark states. Thus in these 



models contradictory results for polarised scattering are obtained (see [11| for a discussion of 
these results) but being in conflict with the OPE their predictions should be ignored. As a 
result we see that there is no remaining theoretical conflict for the case of the covariant quark 
parton model description of polarised scattering. Whether it provides a good phenomenological 
description of polarised scattering is a question for experiment. 



4 Corrections from finite quark masses 

Only for massless quarks is the parton polarisation vector s'^ proportional to k'^. In the m 7^ 
case the partons will, in general, have three polarisation components. Let us choose as a 
convenient basis the vectors rii given by 

< = ^^[(p.A;)P-mV] 

iV2 

iViiV2 

= ^[M\k.S)k^ + NlS^ - {p.k){k.S)p^] (22) 
9 



where 



= [{p.k'f -m^M^]^ 
N2 = -[{p.kf - M''{k.Sf -m^M^]^ (23) 

In the parton rest-frame these vectors form a set of orthogonal unit 3-vectors and, in any 
frame, satisy k.Ui = 0, ni.rij = —Sij. In the proton rest-frame (PRF), they are given by 

ni = -{\k\,kok), ^2 = X 5), ks = ^[k X {k X S)] (24) 
m sm (J sm (J 

where k.S = cos6 and k,S are unit vectors in the direction of the parton momentum and 
proton spin respectively. The charge weighted combination of the relevant parton distributions 
in each direction we call Afi{p, k, s). The vector rii is the longitudinal polarisation vector of 
the parton, mn'^ = k^ + 0{rn? /M'^) and so A/i(p, /c, s) is to be associated with the Af{p, k, s) 
of eq(D. 

These 0{m? / M'^) corrections modify the integrand of eq(§) and since the value of k"^ = 
M'^[x{y — x) — m^/M^] the integrand of eq(|^) is also appropriately modified. In addition, there 
are 0{m'^ / M"^) corrections to the limits of the integration in y. If we write a = m'^/M'^ then 
ymax = l + a follows from constraining {p — kY > while ymin = x + a/x follows from k^ > 0. 
The resulting expressions for gi and (72 associated with the component rii are 

fi)/ Z"^™"'' 2a, ~ , , 
9i [x) = dyx{2x-y){l ) h{y) 

92\x) = dy[x{2y-?>x)+a{A--?>)]hM (25) 

O Jymin y 

Here, Hy) = (1 - ^ dk^ f{k\y) where A/i(p,fc,5) = -{k.S)/M ~f{k\y) 

Notice that ymax > ymm only if x > a, i.e. g^il are zero for x < m? /M"^. The above 
expressions do not satisfy the WW sum rules(|D but g2^ does satisfy the BC sum rule(||) 
exactly. To see this we assume that we can interchange the order of the x, y integrations and 
use the fact that the support of x, for fixed y, is given by Xmin ^ x < Xmax where Xmin and 
Xmax are solutions of x^ — xy + a = 0. 

The component is irrelevent since it is orthogonal to and so Af2{p, k, S) vanishes. 
The rig component does contribute however and substituting into eqs(|^,|) gives the appropriate 
^ Increasing this threshold to positive values of 0{m?) leads to only minor changes since then a; > a + 0{a^) 
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contributions to gi and (?2 

x{y — x) — a 



Jymin 4iVliV2 J 

2 



gf\x) = TrmM' / dy—-— ^ - 3xiy - x) + a dk^AMp^k^s) (26) 



'Vmin SN1N2 

.(3) .(3) 



where ymim Umax and a defined as above. The gi and g2 defined by eq ipq ) satisfy analogous 
sum rules to the WW sum rules in the limit a — > 0. Provided A/3 oc N2, which is true in any 
reasonable model (see below), we can write in this limit 



^ dxx""-^ 



(27) 



An 

This relation is interesting in that not only is the BC sum rule satisfied but also the first 
moment of 5^2^"* vanishes ^ which implies that g^^ is suppressed in general. When a 7^ the 
n = 1, 2 sum rules both survive exactly. 

The above expressions, eq(|25D and eq(^), are quite general. As they involve two distribution 
functions, they do not lead to a relation between gi and g2- However, as we now discuss, the 
covariant parton model allows us to develop a plausible model to describe the source of the 
parton polarisation in the proton and hence to relate Afi{p, k, s) and Af^lp, k, s). How are 
the partons polarised? Ultimately it is due to the interaction with the external magnetic field 
which causes the proton spin, through the interaction of the proton magnetic moment, to align 
with the field. At the parton level, while they too will interact directly with the magnetic 
field, they will also have spin-spin interations which align the individual parton spins. The 
latter interaction, being due to the strong force, may be expected to dominate the relative 
parton spin alignment which is what relates Afi{p, k, s) and Af^lp, k, s). Since, on average, 
any combination of parton spins is proportional to the proton spin we are led to a model in 
which the relative magnitude of A/j is due to the spin-spin interaction between the partons 
themselves, the interaction being proportional to the scalar product of the individual parton 
spin with the proton spin. Thus we have 

A/,(p, k, S) = in,.S) m /{p.k, k') (28) 



From eq(22) we have 

{p.k){k.S) 



m (rii.S) 



Ht has been argued that the same is true too for the contribution involving ghions ||T^, p^ . 
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m {n2.S) = 

N2 

m {n^.S) = -m — (29) 

We see that factoring out [k.S) in eq(^) in the massless case is consistent with this model. The 
factor m on the rhs of eq(^) ensures the desired hmit for the Ui contribution as m — > 0. We 
can now write A/3 in terms of /, 



and substitute into eqs (|26|) to get 



Vmax 



gi^\x) = nrn^ I dy 



x{y — x) — a 



y2 



y2 

^ — 3x(?/ — x) + a 



hi{y) 

hi{y) (31) 



2y' 

Thus gi{x) = gi^\x) + gi^\x) and g2{x) = g2'\x) + g2^\x) are determined by a single function 
once again, as in the m = case. The resultant prediction for g2{x) in terms of the measured 
gi{x) is pursued in the next section. 

In leading order the final form of the quark mass effects consists of two pieces. The first 

2 

which starts at 0{^) comes simply from imposing parton kinematics on the scattering process. 
The second, of O(^), is given in eq(PUD and comes from the transverse polarisation states. The 
question immediately arises whether it is the constituent or current mass that is relevant. At the 
level of the quark parton model itself this cannot be answered as it is the QCD interactions that 
cause masses to "run". However one may make an educated guess as to the most appropriate 
choice. The calculation of the parton model process requires, for gauge invariance, that the 
mass of the intermediate parton be the same as that of the initial parton. The former should 
certainly be taken as the current-quark mass at the scale since the large momentum involved 
in the deep inelastic scattering process flows through it. Hence the struck quark mass should 
also be taken to be the current-quark mass. As a result we expect the corrections of 0{jj) 
in eq(p6D to relate to the current quark mass. The kinematic corrections however have two 
origins. The first, due to the ymin cut and from the form of the parton momentum k, comes 
from putting the struck parton on mass shell and again for consistency should be taken to be 
the current-quark mass. The second, the ymax cut, comes from imposing parton kinematics on 
the intermediate states in the lower bob of Fig lb. This must surely be the constituent mass as 
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no large momenta are involved. However, in practice, the former effects are the most important 
and hence we expect the dominant mass effects to be assocated with the current quark mass (in 
the phenomenological analysis given below we make no distinction between the various masses 
- consistent with the parton model interpretation - but, following this discussion, we expect the 
masses should be of current quark magnitude). 



5 Phenomenological analysis 




Figure 2: Fit to the data on gi{x) from refs(4,5,6) using eq.(6) and the comparison of the resulting 
prediction for g2{x) from eq.(7) with the preliminary data of ref(l). This is the m = case. 

Fig. 2 shows the present situation for the experimental measurements by SLAC[^], EMC]^ 
and SMC0 for gi{x) together with the recent preliminary measurements of g2{x) at SLAC||1|. 
The range of over which all the measurements are taken is fairly wide and, in principle, 
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one should try to take account of possibly sizeable variation at fixed x values. However we 
do not tackle this here as we are concerned here only with the challenge of trying to predict 
the size and shape of (72(2;) from that of g\[x) at some canonical value of where the parton 
model may apply and, by definition, QCD corrections are assumed to be small. 

Also shown in Fig. 2 is the expected m = prediction for g^ix) where we fit to the data on 
g\[x) using the expression eq(^ with an assumed form for h{%j) = ax~^{l — xY which is then 
inserted into eq(|^) to give g2{x). That is, the curve on Fig. 2 is the prediction for (72 (x) given 
by the WW relation eq(0). We find 6 ~ 3 and c ~ 1. Certainly at large x it is encouraging 
to see the preliminary g2 data lying close to the WW expectation. Note that, from eq(|ip, 
the prediction for (72 (x) at any x value depends only on values of x at and above that value 
and so is not sensitive to uncertainties in the small x region. Likewise, parametrisations of 
the distributions h{y) are not required to satisy various theoretical conjectures for the small x 
behaviour since this region is quite irrelevant to our concerns here. 

Our conclusion following from Fig. 2 is that the preliminary (72 data, albeit with relatively 
large errors, is consistent with WW relation and hence, within the framework of the covariant 
partom model, consistent with a light quark mass of zero. The next step is to ask if the data 
are also consistent with a sizeable quark mass. 

To answer this question, we carried out fits to gi{x) = gi^\x) + gi^\x) using eqs. (p5| , |3l|) 
for m/M up to 0.2 and a similar parametrisation for h{y) as in the m = case above. The 
resulting fits, together with the corresponding predictions for (72 (x) = g2^\x) + g^2\^) shown 
in Fig.3. The quality of the fits is good provided m/M < 0.04, however the resulting never 
improves on the value of 59 for 50 data points achieved by the m = fit. This remains true 
even when more complicated parametrisations of h{y) are considered. 

The components of gi{x) and (72(2;) associated with the nl and components of the parton 
polarisation are shown in Fig. 4. As m/M increases, gl (x) grows at low x — the suppression 
of 0{m'^ / M"^) is largely offset by the denominator y"^ in the integrand — and tends to spoil the 
quality of the fits. Notice in Fig. 4 that each component of (72(2;) correctly integrates to zero 
and note how g2'\x) is suppressed relative to g2^^{x) due to the vanishing of the first moment. 

To quantify the violation of the WW sum rules in this particular model for the quark mass 
effects we consider the ratio r„ given by 

^Idxx^-^ [r^g,{x) + g2{x)] 



J^dxx--^ l^giix)] 
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Figure 3: Fits to the data on gi{x) from refs(4,5,6) using eqs. (25,31) for m/M up to 0.2 and the 
comparison of the resulting predictions for g2{x) from eqs. (25,31) with the preUminary data of ref(l). 
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Figure 4: The two components gi^\x) and gl'^'{x) for the fits shown in Fig.3 
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That is, r„ is the WW moment normahsed by the gi contribution to that moment. Up to values 
of m/M where the model is able to successfully describe the data, can be used as a direct 
measure of the quark mass. For values of m/M = 0.01, 0.02, 0.04 we get r2 = 1.2, 3.4, 9.0%. 

As the precision of the SLAC g2 measurement increases, we expect that a phcnomenological 
analysis such as this can offer a new and practical procedure for testing whether the quark 
model is a good approximation. 
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